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ABSTRACT 

We present a model for the Type la supernova remnant (SNR) of SN 1604, also known as Kepler's SNR. We find that its main features 
can be explained by a progenitor model of a symbiotic binary consisting of a white dwarf and an AGB donor star with an initial mass 
of 4 - 5 Mq. The slow, nitrogen rich wind emanating from the donor star has partially been accreted by the white dwarf, but has also 
created a circumstellar bubble. Based on observational evidence, we assume that the system moves with a velocity of 250 kms -1 . 
Due to the systemic motion the interaction between the wind and the interstellar medium has resulted in the formation of a bow 
shock, which can explain the presence of a one-sided, nitrogen rich shell. We present two-dimensional hydrodynamical simulations 
of both the shell formation and the SNR evolution. The SNR simulations show good agreement with the observed kinematic and 
morphological properties of Kepler's SNR. Specifically, the model reproduces the observed expansion parameters (m = V/(R/t)) of 
m « 0.35 in the north and m « 0.6 in the south of Kepler's SNR. We discuss the variations among our hydrodynamical simulations 
in light of the observations, and show that part of the blast wave may have traversed through the one-sided shell completely. The 
simulations suggest a distance to Kepler's SNR of 6 kpc, or otherwise require that SN 1604 was a sub-energetic Type la explosion. 
Finally, we discuss the possible implications of our model for Type la supernovae and their remnants in general. 
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1. Introduction 

Type la supernovae (SNe la) are of prime interest for many areas 
in astrophysics. They are important cosmological standard can- 
dles, because of their high, and well calibrated ( Phillips et al. 
1992j > peak luminosities. This has led to the realization that 



we appear to live in a Universe whose expansion is accelerat- 
ing (Perlmutter et al. 1998; Garnavich et al. 1998 ). In addition, 
SNe la are major contributors to the chemical enrichment of the 
Universe, as they are the principal source of iron peak elements. 

There is a consensus that SNe la are the result of thermonu- 
clear explosions of CO white dwarfs (WDs) in binary systems 
that approach the Chandrasekhar mass due to either accretion 
from a companion star (single degenerate scenario), or by the 
merging of two WDs (double degenerate scenario) (see the re- 
views byjHillebrandt & Niemeyer|2000|[Liviop ~ 



degenerate (DD) s cenerario is favored by stellar po pulation syn- 
thesis models (e.g.lYungelson & Livio|2000[ jRuiter et al."|2 009 ; 



2000). The double 



|Claeys et al.||201l| , but the explosion mechanism itself is far 
from clear. The main problem is that WD mergers lead to an 
off-center ignition, which converts carbon and oxygen into oxy- 
gen, neon and magnesium. This results in an accretion induced 
collapse and the formation of a neu tron star, rather than a ther- 
monuclear explosion (e.g. [Saio & Nomoto||l 985; Mochkovitch 
|& Livio|199"0l|Nomoto & Kondo|199lT ~ 

For the single degenerate (SD) scenario the most puzzling 
problem is the nature of the progenitor binary systems. The WD, 
in order to reach the Chandrasekhar limit, should accrete and 
burn the mat erial from its c ompanion star at a rate of around 
10" 7 M yr _1 (Nomoto 1982). For lower accretion rates the accu- 



mulated material is liable to unstable burning, giving rise to nova 
explosions that expel more mass than is accreted. On the other 
hand, high accretion rates lead to the expansion of the WD's 
photosphere to red giant's dimensions and, in the absence of sta- 
bilizing effects of strong accretion winds (Hachi su~et al.fl 996, 
1999}, the system will undergo a common envelope episode, 
which does not lead to a SD SNe la pben & Tutukov|1984| ). The 
finetuning needed to create a SN la explosion is in sharp contrast 
with the relatively high probability of ~ 15% for stars in the 3-8 
M range to eventually explode as S N la (e.g. Mann ucci et al. 
[20n6HdePlaa et al.|2007HMaoz|2008) . 

Our understanding of SN la would benefit greatly if we 
would identify progenitor systems directly or indirectly through 
the imprints they have on the SN spectra, or the supernova rem- 
nants. Several studies have been conducted along these lines, but 
reach different conclusions about the origin of the SNe la (e.g. 



Ha muy et aLl2003l|Panagia et al.|2006||Patat et al.|2007||Hayden| 
et al.|2010| ). 



Here, we report our investigation of what could be the origin 
of the remarkable structure of the circumstellar medium (CSM) 
that shapes Kepler's SNR (Kepler for short), the remnant of 



the h istorical SN 1604 ( [Green & Stephenson||2003[ |Blair eraL 
2007). As we will discuss below, SN 1604 was very likely a SN 
la and the CSM observed in the evolved SNR puts constraints on 
the type of progenitor of this SNR. 

Kepler (G4.5+6.8) is located relatively high above the 
Galactic plane, 590^5 pc, with d$ the distance in units of 5 kpc. 
Its radius is 2.6ds pc. The distance itself is not well known. 



* email: a.chiotellis@astro-uu.nl 



Reynoso & Goss (1999 ) used the HI absorption feature to place 
a lower limit of (4.8 ± 1.4) kpc and independently placed an up- 
per limit of 6.4 kpc based on the lack of absorption by an HI 
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250 kms -1 (Blair et al 



Fig. 1. Chandra X-ray image of Kepler's SNR, with red indicat- 
ing Si-Ka emission (1.75-1.95 keV), green Fe-L emission (0.8- 
1.6 keV), and blue continuum emission (4-6 keV). The image is 
based on a deep, 750 ks, Chandra observation ( [Reynolds et al. 
12007] ). 



cloud. Sankr it et al.| ( |2005| ) combined an estimate of the shock 
speed based on the Ha line width with a proper motion mea- 
surement to derive a distance of d = 3.9^ kpc. Although these 
distance measurements agree with each other within the errors, 
the lack of a detection of Kepler in TeV gamma-rays by H.E.S.S. 
( Aharonian et al. 2008 ), coupled with gamma-ray model pre- 
dictions ( [Berezhko et al.|[2006| ) and the energet ics of the SNR 
based on expansion measurements (Vink 2008 ), suggest a dis- 
tance > 6 kpc, or otherwise a subenergetic explosion. 

The SNR shows bright optical nebulosity with prominent 
[Nil] line emission in the north, indicating dense material with 
elevated nitrogen abundan ces, [N]/[N] Q > 2, b ut otherwise 
roughly solar metallicity ( [Blair et al.|1991j ) . Radio (|Dickel et al.| 
1988) and X-ray expansion measurements dVink 2008; Katsuda 
et al. |2008| ) of the SNR show an overall expansion parameter 
m = V/(R/t) « 0.6, with V the plasma or shock velocity, R the 
corresponding radius, and t the age of the SNR. An exception is 
the northern region, where the expansion parameter is m « 0.35, 
lower than expected for the Sedov-Taylor phase in a homoge- 
neous medium (m = 0.4). Based on these results, Vink ( |2008| ) 
estimates that the nit rogen-rich shell must have had a mass of at 
least 1 M Q (see also Bla ir et al.|2007) . 

The presence of the nitrogen-rich shell has puzzled as- 
tronomers for l ong and has led to the claim that SN 1604 was 
a Type lb SN ( |Bandiera|[T987l ), with the shell being the shed- 
ded outer envelope of the progenitor. In order to explain the 
height above the Galactic plane, the one-sided morphology of 
the shell and the high proper motion of the nitrogen-rich knots, 
Bandiera ([1987 ) argued for a high proper motion of the progen 



itor, ~ 280 kms (see also Borkowski et al.||1992). This high 



systemic velocity of Kepler's progenitor was verified observa- 
tionally based on the proper motion and radial velocities of the 
nitrogen-rich knots of the remnant (Bandiera & van den Bergh 
|1991| ) and by the H a narrow component of the non-radiative 



shocks of the SNR, which gave u* 
[T99T]|Sollerman et aL|2003) 

Since the late nineties it has become clear that SN 1604 was 
very likely not a Type lb, but a SN la. The main reason is that 
the X-ray spectrum has prominent Fe-L emission and relatively 
little oxygen emission (Kinuga sa & Tsunemi||1999[ [R eynolds 
|etal.|2007| ), characteristic of Type la SNRs ( [Hughes et al.|1995] h 
Supporting evidence for a Type la identification is the presence 
of B aimer dominated shocks and the absence of an X-ray emit- 
ting, cooling, neutron star ( [Reynolds et al. [2007 ). This means 
that the dense circumstellar nitrogen-rich shell must originate 
from the progenitor system of the SN la. 

Here, we show that the characteristics of Kepler can best be 
explained within the framework of a single degenerate SN la 
model, assuming non-conservative mass transfer through wind 
accretion from a 4-5 M asymptotic giant branch (AGB) star. 
We retain in this model the idea that the northern shell is the 
result of a bow shock caused by the motion of the progenitor 
system (Bandiera 1987; Borkowski et al. 1992), and we adopt 
the observed systemic velocity of 250 km s" 1 . Hydrodynamical 
simulations show that this model can account for the morphol- 
ogy of the SNR and its expansion characteristics. 

In Sect. [2] we discuss the progenitor binary system and in 
Sect. [3] we determine its implications on the requirements for 
the interstellar medium (ISM) and the evolution of the circum- 
stellar medium and SNR. In Sect. [4] we model the system using 
hydrodynamic simulations with the appropriate parameters and 
discuss the differences when varying the exact parameters. We 
evaluate the results of the simulations and our progenitor model 
in a broader context in Sect. [5] and end with our conclusions in 
Sect. 



2. A type la progenitor scenario for SN 1604 

A model of Kepler's progenitor system should explain the for- 
mation of a > 1 M Q asymmetric shell with solar metallicity and 
enhanced nitrogen abundances, which lies at the northern region 
of the remnant at a distance of 2-3 pc from the explosion center. 

Potentially, this shell could have been formed by i): substan- 
tial outflows of the WD's surface such as nova explosions or 'ac- 
cretion winds', ii)\ the wind of the WD progenitor star (i.e. the 
SNR interacts with a planetary nebula-like shell), Hi): the wind 
of the donor star, or by iv): the ejected common envelope, in the 
case of a DD progenitor scenario. 

The nova explosions, which are related with SNe la are the 
reccurent novae that occur on the surface of massive WDs ( 
> 1.2M G ) in the last phase of the binary evolution (Hachisu 



& Kato]|2001| 



Hachi su et al. 2008| ). However, during the rec- 



curent novae phase, the total mass that is ejected from the WD 
surface is of the order of 10 -3 - 10- 2 M o ( [Hachisu etal.|2008 ). 
Giv en that recurrent nova ej ecta have similar to solar abundances 
1992[ for summary), they are not able to 



(see|Livio & Truran 



accumulate enough heavy elements into the CSM shell to re- 
produce the observed chemical abundances. On the other hand, 
an outflow in the form of an accretion wind emanating from 
the WD surface is so fast (u win d ~ 1000 kms" 1 ) that it would 
form a large low-den sity cavity around the progenitor system. 
Badenes et al. (2007) showed that these cavities are at odds with 
the observed radii, shock velocities and ionization timescales of 
Galactic, LMC and M31 Type la SNRs (among them Kepler 
SNR was included). Our hydrodynamical simulation of a SNR 
evolution in a cavity formed by accretion winds verifies their re- 
sults. 
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0.6 
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Table 1. The production factor and the chemical abundances as 
a function of the stellar (initial) mass o f an AGB star with solar 
metallicity (Karakas & Lattanzio 2007 ). For the case of 5 M© the 
AGB mass loss rates by Vassiliad is & Wood ( 1993 ) have been 
used. 



Case ii can be excluded because of the time scales involved. 
The time interval between the planetary nebula phase and the 
Type la exposion should be at least ~ 10 6 yr as the WD needs to 
accrete at least 0.2 M© at a rate of 10" 7 M yr -1 . In such a period 
the formed shell would have collapsed under the ram pressure of 
the interstellar medium as the binary system is moving with a 
velocity of 250 kms" 1 . The same line of argumentation can be 
applied for case iv, since, after the ejection of the common en- 
velope the merging timescale of the two WDs is at least 0. 1 Myr 
(Claeys 2011, private communication). Although this scenario 
seems unlikely, the lack of understanding of the ejec tion of a 
common envelope and the subsequent evolution (e.g. |Taam &| 
Ricker 2010) prevents us from drawing a definitive conclusion 
about the likelihood of case iv as a progenitor scenario. 

Therefore, the wind from the donor star (case Hi) appears to 
be the most likely origin for this circumstellar shell as the wind 
velocities of evolved stars are much smaller than those of the 
WD's outflows, resulting in smaller and denser cavities and their 
formation sustains until the moment of the explosion. 

Chemical composition of the circumstellar shell 

The northern shell of Kepler reveals a roughly solar metallic- 
ity plasma with elevated nitrogen abundances of at least 2 times 
the solar value. The only process to enrich the surface of a low 
mass star with nitrogen is hot bottom burning (HBB), which oc- 
curs in massive AGB star s (M > 4 M©, depending on metallicity, 
Karakas & Lattanzio 2007). During this process the convective 
envelope of the star dips into the outer layer of the H-shell, which 
results in a thin layer hot enough to sustain proton-capture nucle- 
osynthesis. HBB converts 12 C into 14 N and if the third dredge-up 
is operating, the stellar atmosphere will become nitrogen rich. 
Based on the results of Karaka s~& Lattanzio] ( |2007| ), the donor 
star of the progenitor system was most likely an AGB star of 
4 - 5M with solar metallicity (see Table [TJ. 



Symbiotic binary 

Efficient mass accumulation on the WD requires stable mass 
transfer. Roche-lobe overflow (RLOF) by stars with deep con- 
vective envelopes, such as the AGB donor star of our scenario, 
leads to dynamically unstable mass transfer and the formation 
of a double WD binary system engulfed by their common enve- 
lope. Such a system is not possible to end up as a SD Type la 
SN. 

Thus, the progenitor system requires a wide symbiotic bi- 
nary, where, through non-conservative mass transfer, part of the 
donor's stellar wind is accreted onto the WD, while the rest is 
ejected from the system forming a circumstellar shell. Since 
a CO WD needs to accrete at least 0.2 M in order to ap- 
proach the Chandrasekhar mass, and the envelope mass of the 
AGB star is ~ 3-4 M©, an (average) accretion efficiency of 
> 5 — 10 % is required for a Type la explosion. The mass trans- 
fer can progress by either pure wind accretion, or by the more 



promising and e fficient method of wind-RL QF accretion as sug- 
gested by Mohamed & Podsiadlowski ( 2007 ). In either case, the 
AGB star should stay within its Roche lobe, implying that the 
orbital radius has to be larger than 1600 R© for typical mass ra- 
tios. 

Systemic motion 

Considering that the progenitor binary was moving with a 
systemic velocity of 250 kms -1 , Kepler' s Galactic latitude 
and the asymmetric accumulation of the circumstellar shell can 
be readily explained by its supersonic motion away from the 
Galactic plane (Bandiera |1987[ Borkowski et al.|1992 ). The in- 
teraction of the wind with the ram pressure of the interstellar 
medium causes the formation of a bow- shaped shell. The subse- 
quent supernova drives a blast wave that currently interacts only 
with the nearest part of this shell. 

3. Formation and evolution of the wind bubble and 
of the subsequent SNR 

Now that we have outlined the properties of the progenitor 
binary system, we can study the formation and evolution of the 
wind bubble and the subsequent interaction of the SNR with 
it. In this section, we approach the dynamics analytically and 
specify the variables, which determine our model. 



3.1. The formation of the wind bubble and the bow shock 
model 

The properties of a supersonically moving wind bubble, interact- 
ing with the ISM, is a function of four variables; the mass loss 
rate, M, the wind velocity u w , the systemic velocity, w*, and the 
ISM density, n{ sm . The system will be in a steady state, when the 
timescale of the interaction process between the wind and the 
ISM will be larger than the flow timescale defined as: 



rflOW "(i)^ 10 li)(l^^) lyr ' 



(1) 



where r in this case indicates the distance of the bow shock from 
the mass-losing star ( [Borkowski et al.|1992| ). 

In a steady state, the position of the bow-shock shell is de- 
termined by the balance of the stellar wind and the ISM ram 
pressure. At the so-called stagnation point the wind termina- 
tion shock has the shortest distance, ro, to the mass losing star. 
Equating the momentum fluxes of the ISM and the CSM, this 
can be estimated to be around: 



r = 1.78 x 10 3 



M w u, 



1/2 
2 I PC, 



(2) 



when the vel ocities are measured in km s" 1 , M w in M yr _1 and 
n ism in cm" 3 ([Huang & Weigert|l982l ). 

In Sect. 2|we argued that the donor star of the progenitor sys- 
tem was a 4 - 5 M© AGB star. These stars are characterized by 
dense, slow, radiation-driven stellar winds with typical terminal 
velocities of 5 - 20 km s" 1 and temperatures of 1000 - 2000 K 
(Habing & Olofsson 2003 ). The mass loss rates are fluctuating 
in a wide range, starting from 10" 8 - 10" 6 M yr _1 during 
the early AGB phase to 10~ 6 - 10" 5 Mpyr" 1 at the thermal 
pulsating AGB phase (Vassiliadis & Wood 1993[ |Karakas & 
[Lattanzio 2007 ). The total duration of the AGB phase is on 
the order of 1 Myr. Assuming typical wind loss parameters 
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Fig. 2. The mass loss rates versus the wind velocity of the AGB 
star for the values which satisfy r = 2 pc and r = 3 pc for 
fticm = 0.5 x 10" 3 cm -3 and u* = 250 km s" 1 . 



appropriate for the thermal pulsating AGB phase, a steady state 
situation is reached after ^ ow « 0.1 - 0.3 Myr for ro = 2 pc and 
^flow ~ 0.1 - 0.4 Myr for r = 3 pc (Eq.[T]), which is well within 
the limits of the lifetime of an AGB star. 

In order to reproduce the radius of the stagnation point, we 
also need an estimate of the ISM density. We do not have a pri- 
ori estimates of this density, as there are no clear observational 
data for its properties in Kepler's neighborhood. Both the neu- 
tral/ionized warm and the hot ionized component extend to its 
latitude (Reynolds 1991). However, the constraint on the dis- 
tance of the stagnation point (r s 2 - 3 pc, see eq. [2]) requires 
rii sm < 10~ 3 cm~ 3 . This value is characteristic for the hot ion- 



ized component (T ~ 10° K), where 



10-4 _ lO- 3 cm" 3 



(McKee & Ostri ker|[i977| ) and was also used by Borkowski 
|et al.| ( |1992| ), basedo n the observed densities of the win d shell 
in she u « 1 - 10 cm- 3 , |Smith et al.|l989[[Blair et al-|l99l} - 

Figure[2]shows the possible combinations of M w and u w that 
result in ro = 2 pc for the lower distance estimate of Kepler, and 
r = 3 pc for the larger distance estimate. In this plot we have 
adopted n{ sm = 0.5 x 10" 3 cm -3 for the hot ionized ISM, and 
u* - 250 km s" 1 . 



3.2. Supernova properties and evolution 

The canonical values for the mass and energy of Type la SNe are 
Af e j = 1.4 M© and £ e j = 10 51 erg (Woosley et al.|2007|), resp ec- 
tively. We apply the self- similar solution ol Chevaliefp f 1 982| ) for 
the case where the SNR expands into an p oc r~ L wind profile, in 
our model corresponding to the regions interior to the shell. The 
ejecta is assumed to consist of a constant density core with an en- 
velope that follows p oc r -7 density profile, while the velocity in 
the ejecta increases linearly. By assuming energy conservation, 
the expansion of the ejecta in a wind bubble with p w ind = qr~ s 
with s = 2 and q = M w /(4nu w ) is given by: 



13x[Ag n /q] lKn - s) t (n - 3)/(n - s \ 



(3) 



where A is a constant equal to 0.27, and t is the age of the rem- 
nant (in sec). Finally g is a constant given by: 

g 1 = (25/21;r)(4/M ej ). 

For the age of Kepler's SNR of t & 400 yr, mass-loss param- 
eters of M w = 10" 6 - lO^Moyr" 1 , and u w = 10 - 20 km s"\ 



we find R snr = 3.0 - 4.7 pc. These values of the SNR radius 
correspond to a distance of Kepler around d = 6.2-9.7 kpc. 
In order to obtain a SNR radius consistent with a distance of 
d = 4 kpc, we also have t o consider a subenergetic explosion 
of E = 0.2 x 10 51 erg (c.f. |Vink||2008] ). In that case we obtain 
R sm = 1.6 - 2.4 pc and d = 3.3 - 4.9 kpc. 



4. Hydrodynamic modeling 

We employ the hydrodynamics code of the AMRVAC frame- 
work ( [Keppens et al.|2003] ) for our simulations of the circumstel- 
lar bubble around Kepler's progenitor system and the subsequent 
evolution of the supernova ejecta. We perform the calculations 
on a 2D grid in spherical coordinates and assume symmetry in 
the third dimension. The Euler equations are solved conserva- 
tively with a TVDLF scheme, using the adaptive mesh strategy 
to refine the grid where needed as a result of large gradients in 
density and / or energy. Our radial span is 2 x 10 19 cm and the 
range of the polar angle is 0° to 180°. On the base level, we use 
240 x 120 cells (R x 6) and we allow for three refinement levels 
during wind evolution and four for the SNR evolution, at each of 
which the resolution is doubled. The maximum effective resolu- 
tion, thus, becomes 1.04 x 10 16 cm by 0.188°. Radiative cooling 
is prescribed using the cooling curve by Schur e^etaLl ([2009). 

First we simulate the formation of the CSM bubble with the 
ISM bow shock shell. We model the system in the rest frame 
of the progenitor system and model the ISM interaction as an 
inflow. The ISM with density pi enters the grid antiparallel with 
the y-axis and with a momentum m r = p[U* cos 6 (see Fig. [3}. 
Thus, the symmetry axis is aligned with the systemic direction of 
motion, roughly corresponding to the northern region of Kepler's 
SNR. In the inner radial boundary, we impose an inflow in the 
form of a cold, slow, stellar wind with a density profile of p = 
M w /(4nr 2 u w ) and momentum m r = pu w and m e = 0. In the 
second stage, we introduce the supernova ejecta into the wind 
bubble and allow the SNR to evolve. 

In Sect. [3] we defined the constraints for the systemic veloc- 
ity, mass-loss rate and wind velocity. These constraints allow for 
some variation in those parameters, which will all give roughly 
similar values for the radius of the shell, but which may never- 
theless influence the properties of the SNR. In order to inves- 
tigate this, we did several hydrodynamical simulations, varying 
the wind parameters and systemic velocity. In addition, as there 
is some uncertainty concerning the distance to Kepler, we study 
the cases, where the stagnation point is placed at ~ 2 pc and 
at ~ 3 pc from the explosion center, corresponding to the SNR 
distances of 4 kpc and 6 kpc respectively. The grid of models is 
summarized in Table [2] 



4.1. Wind bubble evolution 

Figure |4]illustrates the evolution of the pre-supernova wind bub- 
ble using model A (see Table [2|. It shows the typical four-zone 
structure of a stellar wind bubble, consisting, from inside out, 
of the freely streaming wind with p oc r~ 2 , the shell of shocked 
stellar wind, the shell of shocked ISM, and the unperturbed in- 
terstellar gas. At the border of each 'zone' a density jump occurs 
due to the termination shock, the contact discontinuity, and the 
outer shock. The shear flow at the interface of the two fluids is 
susceptible to the Kelvin-Helmholtz instability, which results in 
a wavy structure of the contact discontinuity. 

The wind termination shock initially propagates rapidly out- 
wards, but then decelerates, until it is stabilized, due to the mo- 
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Table 2. The properties of the 5 studied models. Additional parameters that have been used for all the models are: T w = 
1000 K, r ism = 10 6 K and M ei = 1.4 M . 
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Fig. 3. The 2D density profile of the simulation of the bow- 
shocked wind bubble. The mass-losing star is located in 
the origin where we let a radial wind flow enter the grid. 
Simultaneously, isotropic, homogeneous gas is entering the grid 
antiparallel with the y-axis, which represents the motion of the 
ISM in the star's rest frame. The arrows correspond to the vec- 
tors of the momentum of each flow. 



mentum equilibrium between the wind and the ISM. This equi- 
librium is first achieved at the stagnation point, and then propa- 
gates to larger azimuthal angles. 

We introduce the SN explosion once the stagnation point is 
at 2-3 pc, while taking care so that the donor star's mass loss 
is less than the mass of the envelope of the AGB star and the 
duration is less than the life time of an AGB phase. In model A, 
we find that these conditions are met at a time t = 0.38 Myr after 
initializing the evolution of the bubble. The stagnation point in 
this case is located at a distance of ~ 2.7 pc, the outer shock 
is at a distance of 5.2 pc, and the donor's mass loss is equal to 
Mt = 3.8 M . 

4.2. SNR evolution 

Fig. [5] shows the evolution of the SNR for model A in an 
ambient medium that corresponds to the wind bubble at time 
t = 0.38 Myr. Initially, the SNR evolves in the spherically 



symmetric, freely expanding, wind-blown bubble (left panel). 
The freely streaming ejecta, shocked ejecta, contact disconti- 
nuity and the shell of the shocked CSM can be distinguished. 
Rayleigh-Taylor instability develops at the contact discontinuity 
between the ejecta and the shocked CSM. The thin black lines 
outline the interaction regions between the ejecta and the CSM, 
while the white ones indicate the interface between the CSM and 
the ISM. Around ~ 285 yr after the explosion (second panel) the 
blast wave starts to interact with the shocked stellar wind of the 
circumstellar shell in the area of the stagnation point. The blast 
wave sweeps up the dense shell's material and the deceleration in 
this area is therefore stronger, something that breaks the spheri- 
cal symmetry of the SNR. 

With time, a larger portion of the SNR starts to interact with 
the circumstellar shell. In the snapshot that corresponds to the 
age of Kepler's SNR (t = 412 yr, right panel in Fig. [5]) more than 
one third of the remnant interacts with the wind shell. Around 
the stagnation point, the blast wave has slightly penetrated the 
shell. At this moment, the blast wave in the southern region is 
located at a distance of Ruast ~ 3.6 pc, and the contact discon- 
tinuity lies at 7? cd « 2.7 pc. The highest densities are found in 
the area around the stagnation point, where it reaches values of 
n « 4.0 cm -3 . This area is expected to have the highest emissiv- 
ity. 

Fig. [6] shows the density, expansion rate, pressure, and tem- 
perature of the simulated SNR at the age of Kepler. The ex- 
pansion rate, as expected, is much lower in the region, where 
the interaction with the circumstellar shell takes place. Around 
the stagnation, point we find values as low as ~ 0.35, while for 
the non-interacting areas, the expansion rate is 0.6. These val- 
ues correspond to the expansion rate of the plasma and are con- 
sistent with results from X-ray observations ( Vink 2008). The 
third panel shows the pressure, which clearly marks the shocked 
plasma being rather homogeneous in the entire SNR, consis- 
tent with the expectation of pressure equilibrium between the 
shocks. The stronger interaction in the region around the stag- 
nation point causes the SNR to be thinner and to have higher 
pressure compared to the rest of the remnant. The temperature 
plot (right panel) is limited to the range log T = 6.0-9.5, cor- 
responding only to the shock-heated SNR shell. Note that X-ray 
observations usually provide electron temperatures, which may 
be lower than the overall plasma temperatures depicted here due 
to non-equilibration of ions and electrons. 

Fig. [7] shows the position and the velocity of the blast wave 
(upper plot) and its expansion rate (lower plot) for two differ- 
ent azimuthal angles of 6 = 0° and 6 = 180°. During the evo- 
lution of the SNR in the freely-expanding wind region, these 
three quantities are identical for the two angles. After 300 yr 
the collision with the shell starts and at the interacting region 
(0 = 0°) the shock wave is strongly decelerated, while the SNR 
becomes slightly aspherical. The expansion rate is ~ 0.8 for the 
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Fig. 4. The evolution of the wind bubble of model A. The snapshots from left to right correspond to the times 0.10 Myr, 0.29 Myr, 
0.38 Myr and 0.57 Myr. 
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Fig. 5. SNR evolution of model A. The snapshots from left to right correspond to the times 158 yr, 285 yr, 349 yr and 412 yr. 
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Fig. 6. Model A SNR properties at 412 yr. From left to right: density, expansion rate, pressure and temperature. 
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areas, which are still within the unshocked-wind region, while 
for 6 = 0° it drops to ~ 0.5 for the SNR at Kepler's age. 
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Fig. 7. Model A blast wave properties. The upper plot shows the 
position (solid lines) and the velocity (dotted lines) of the shock 
versus the time. The lower plot displays the evolution of the ex- 
pansion rate. The black color corresponds to the azimuthal angle 
of 0° while the red for 1 80° . 



4.3. Comparison of the different models 

In this section, we evaluate differences in the morphology of the 
SNR, when varying the mass loss properties and the systemic 
velocity of the system, as summarized in Table [2] 

Wind bubble 

In Fig. [8] we show the densities of the bubbles for the different 
models, at the moment immediately prior to the SN explosion. 
The left panel shows model A, which was already described in 
section [47TT 

In model B we use a slightly lower mass loss rate for the 
formation of the wind bubble. Due to the longer duration of the 
shell formation, the Kelvin-Helmoltz instability at the contact 
discontinuity has developed more prominently. The chosen wind 
velocity is slightly higher, while the ISM density is lower in or- 
der to get the stagnation point at a distance of ~ 3 pc. The wind 
termination shock and the contact discontinuity are farther from 
the star than in model A. This model has been considered in or- 
der to retain the blast wave of the SNR within the shocked wind 
shell at the current age of Kepler's SNR. 

For model C we use an even lower mass loss rate and velocity 
of the AGB wind. These values fit better with the Reimers model 
for AGB mass loss parameters ( Reimers||l975 ), or an AGB at 



the early phase, according to the Vassilia dis & Wood ( |1993 ) 
description. In order to keep the stagnation point at a distance 
of 3 pc, also the systemic velocity and the interstellar medium 
density have been decreased. We consider a longer duration of 
the shell formation/mass transfer phase in order to accumulate 
enough mass and let the shell evolve to a distance of ~ 3 pc. The 
lower mass loss rates and wind velocities yield a lower momen- 
tum in the wind material, resulting in a weaker interaction and 
thus a thicker and comparitively more tenuous shell. 

The final panel shows the bubble for the model D. The wind 
mass loss rate and duration of the shell formation are similar to 
those of model B. However, a lower wind velocity and a higher 
ISM density is applied. This causes the shell to be closer to the 
binary and denser. This shell model has been considered in order 
to be able to include a case, where Kepler is located at a distance 
of 4 - 5 kpc, for which we need the shell at a radial distance of 
2.0- 2.6 pc. 

SNR 

The location, density, thickness and radius of curvature of the 
shell all determine the morphology of the SNR upon interaction 
with the ejecta. Figure|9] shows the resulting densities and expan- 
sion rates of the SNR at an age of 412 yr, similar to the current 
age of Kepler. 

As we noted in Sect |42j in Model A the SNR blast wave 
reaches slightly beyond the shell in the region around the stag- 
nation point, and accelerates again as it protrudes. The expansion 
rate of the plasma in the interaction region is consistent with the 
observed values and it reaches a minimum (~ 0.33) just behind 
the blast wave. This region is characterized by the highest den- 
sities and thus, the emissivity is expected to be also relatively 
high. 

Apart from the protrusion through the shell, model B yields 
similar properties. At the snapshot, which corresponds to the cur- 
rent age of Kepler, almost one third of the remnant interacts with 
the wind shell. The radius of the SNR is larger (4 pc compared to 
3.8 pc in Model A), and the minimum expansion rate is slightly 
higher due to the lower densities of the wind region and wind 
shell respectively. 

These effects are even stronger for the case of model C. Due 
to the much lower mass loss rates, the radius of the remnant is 
4.1 pc, while the minimum expansion rates tend to the value of 
0.4. However, a similar morphology with the other two models 
has been formed also in this case. 

In Model D the asymmetry is the highest. The interaction 
is strong because of the proximity and higher mass of the wind 
bubble, which was designed to accommodate the small distance 
estimate of 4 kpc to Kepler's SNR. The blast wave reaches well 
beyond the extent of the shell resulting in a high expansion rate 
(> 0.7) in that area. This, together with the high degree of asym- 
metry, makes this model inconsistent with the observations of 
Kepler. 

In all of these models, the radius of the SNR is fairly con- 
sistent with the larger distance estimate (~ 6 kpc) of Kepler's 
system. We find that in order to explain the observed morphol- 
ogy in case Kepler is as close by as 4 kpc, we need a subenergetic 
explosion. We adopt the parameters of model D for the bubble, 
so as to have a wind shell relatively close to the star, while de- 
creasing the explosion energy to one fifth of the canonical value 
in model Dsub. The radius of the SNR then becomes ~ 2 pc, and 
the asymmetry and expansion rates are consistent with Kepler's 
SNR. 

From these models, it is clear that the morphology does not 
depend on fine-tuning of parameters. For a reasonable range of 
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values for all of the parameters, the expansion rate, size, den- 
sity, and morphology of Kepler's SNR are well reproduced. The 
different models result in some variation in the blast wave posi- 
tion and locations of maximum density and minimum expansion 
rates. 



5. Summary and discussion 

5.1. Comparison of Kepler's SNR with the models 

We have argued above that the kinematics, morphology and 
abundances in the northern region of the Type la SNR Kepler 
are best explained by a model, in which the progenitor system 
was a symbiotic binary with a WD accreting mass from the slow 
(10-20 kms" 1 ) wind from an AGB donor star. The accretion 
efficiency is moderate and most likely of the order of 10-50%, 
which characterizes the wind, or the wind-RLOF mass transfer. 
This allows for sufficient wind luminosity, needed for creating 
a shell, while, at the same time, it allows the WD to accrete at 
a sufficient rate to reach the 1.38 M within the lifetime of the 
AGB phase of the donor star. The initial mass of the AGB star 
must have been 4-5 M Q , in order to explain the enhanced nitro- 
gen abundances in the northern part of the SNR. 

The wind emanating from the system creates a bubble, 
which, due to supersonic systemic motion of the system (250 
kms -1 ), is surrounded by a one-sided shell, created by the bow 
shock. This explains both the morphology of Kepler's SNR, and 
its height above the Galactic plane of 590^5 pc. In this sense our 
model is reminiscent of the stellar runaway model for Kepler 
proposed by Bandiera (1987]) and [Borkowsk i et aT] ([T992), but 
applied to a Type la progenitor system. 

The parameters of the model are well constrained, as the dis- 
tance between the progenitor system and the wind shell depends 
on a combination of the systemic velocity, wind velocity, mass- 
loss rate and ISM density (Eq.[2| Fig. [2]). At the same time, the 
mass-loss rate should be high enough to allow for enough accre- 
tion onto the WD within the AGB phase, while the total mass 
that is lost from the system is limited by the envelope mass of 
the donor star, -3-4 M . A final constraint is that the wind 
velocity should be high enough to allow for the creation of the 
shell within the lifetime of the AGB phase (Eq.[T]). 

We tested our model with several hydrodynamic simulations 
of both the creation of the stellar wind bubble and the subsequent 
SNR evolution, varying the systemic velocity, the wind velocity, 
and mass loss rate. Since there is some uncertainty about the 
distance to Kepler's SNR, and hence about the size of the wind 
shell, we included also two models (model D and Dsub ) result- 
ing in a stagnation point at 2 pc instead of 3 pc. 

All models reproduce the overall characteristics of Kepler's 
SNR; they have one-sided shells, and the expansion parameter of 
the SNR inside this dens e shell drops to m = 0.3 - 0.4, in agree- 
ment with observations ( |Dickel et al.|1988t[Vmk|2008t[Katsuda| 
et al. 2008] ). An exception is model D. This model has the stagna- 



tion point of the wind at 2 pc, relevant for a distance of Kepler's 
SNR of ~ 4 kpc. The shock of the SNR has moved completely 
through the dense shell after 400 yr, and has an expansion param- 
eter of m > 0.7, inconsistent with the observations. In this case, 
in order to match the model with the observations, the SN explo- 
sion energy has to be E = 0.2 x 10 51 erg (model Dsub), which 
was also suggested by Vink (2008 ). Since normal SNe la have 



models therefore favor the larger distance estimate of Kepler's 
SNR, d « 6 kpc. 

Although the overall morphology and kinematics of Kepler's 
SNR can be reproduced, there are some interesting differences 
between the models. Most notable is the fact that in model A 
the blast wave has just emerged out of the shell near the stagna- 
tion point. The question is now whether this is true for Kepler's 
SNR as well. This question is related to the orientation axis of 
Kepler's SNR with respect to that of our model. Examining the 
available literature and X-ray data we think that a good case can 
be made that the north is the projected direction of systemic mo- 
tion of the progenitor system, and that the shock wave in Kepler 
has indeed just emerged out of the wind-created shell. 

The first reason is that the nitrogen rich knots are found 
mostly in the northwest, and also in the northeast, but in the north 
they are almost absent (Fig. 1 and 2 in Blair et aE]|1991| ). Non- 
radiative shocks associated wit h fast shocks in lower density gas 
are more prevalent due north (Sankrit et al. 2005). Similarly, 
more 4-6 keV X-ray continuum emission associated with the 
shock front (blue in Fig. [T]) appears to be present in the north 
than in the northeast and northwest. Most of the X-ray con- 
tinuum emiss ion in Kepler's SNR s eems to be synchrotron ra- 
diation (Reyn olds et"aL][2007] [Vink 2008 ), which requires fast 
shock velocities V x > 2500 kms 'M Aharonian & Atoyan||l999 



Zirakash vili & Aharon ian 2007). This suggests that the shock 
has picked up speed in the northern region. 

Additionally, the H a images show a much thinner region of 
the emission due north compared both to the rest of the 'shell' 
and with the X-ray images. This could well be due to the fact 
that in the ISM hydrogen is completely ionized and no B aimer 
shocks can be expected beyond the shell. This would be an ad- 
ditional argument in favor of the shock having penetrated the 
shell. However,in order to make a firm conclusion based on this 
feature a more detailed comparison between the X-ray and H a 
images is required. Although, in our view, model A seems to re- 
produce Kepler's SNR best, slight adjustments of the parameters 
of model B and C may also result in a blast wave extending just 
beyond the shell. 

Finally, in our comparison of the models with Kepler's SNR 
we would like to point out that the SNR shell in the models in 
the north is thinner than in the south. This brings also the contact 
discontinuity of ejecta and swept up matter closer to the shock 
front. This, indeed, seems to be the case in Kepler's SNR be- 
cause the emissivity profiles corresponding to different ejecta el- 
ements show that in the northern region there is less stratification 
of the ejecta and the eje cta are closer to the forward shock than 
in the southern region ( Cassa m~Chenai' et~aL]|2004| ). A similar 
situation seems to arise in Tycho's SNR, where the proximity 
of the ejecta to the shock front has been attributed to a change 
in the effective equation of state as a result of efficient cosmic 
ray acceleration ( W arren et al.|2005| ). Our simulations, however, 
can explain this, at l east partially, based on pure hydrodynamics 
alone (c.f |Kosenko et al.|2010| ). 

5.2. Implication of the model for the progenitor of SN1604 

The symbiotic binary model that we advocate here for Kepler 
has some interesting consequences for the current state of the 
donor star. So far, the identification of a donor star of a histori- 



explosion energies in excess of 10 51 erg ([Woosley et al. 



Baade 



and given the high iron content and historical light curve 
|1943| ), SN 1604 must have been a fairly normal Type la. Our 



2007| 



cal supernova h as only been claimed for Tycho/SN 1572 ( [Ruiz- 
Lapuente 2004 ), which was mainly based on its large proper mo- 
tion. However, this result is still controversial, since the large 
proper motion of the candidate star is not accompanied by a 
high spin velocity as would be expected for close Roche-lobe 
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overflowing binary progenitors (Kerzendor f et al.|2009| ). In our 
model no such high orbital and spin velocities are expected for 
the donor star because we require a wide symbiotic binary. It 
should nonetheless maintain the high systemic velocity of the 
progenitor system of ~ 250 km s" 1 , something that constitutes a 
clear signature for its identification. 

The donor star of Kepler is expected to be, according to our 
model, an AGB star with initial mass of 4-5 M©, which has lost 
almost its entire envelope due to the mass transfer and the sub- 
sequent coll ision with the blast wave of the supernova ( Marietta 
et al. 2000). After the collision and the mass strip effect, the 
donor star will expand to re-establish its thermal equilibrium. 
Stars with convective envelopes are characterized by short ther- 
mal timescales, which means that the donor star has most likely 
already re-expanded to its original size. However, further investi- 
gation is needed before a firm conclusion can be be drawn about 
the present-day characteristics of the donor star, as the expan- 
sion timescale is determined by the mass of the stripped layers 
and the energy deposited at the rema ining layers of the do nor 
star due to the collision with the SNR (Podsiadlowski 2003 ). 

Three possible cases can be distinguised: i) The donor star 
has re-established its thermal equlibrium: In this case the donor 
star should be an evolved star with a bolometric luminosity of 
an AGB star with initial mass of 4-5 M Q . The reason for this is 
that giants' luminosity is mainly determined by the mass of the 
core, which remained unaffected by the collision. We, therefore, 
expect an absolute magnitude of M v = -4.5 + 0.5, wh ich at a 
distance of 6 kpc and Ay = 2.8 (Reynolds et al. 2007 ) implies 
an apparent magnitude of my = 12.0 + 0.5. Due to the mass loss, 
the surface gravity of the donor should be much lower that this 
of a 4-5 M AGB star and perhaps the remaining envelope has 
picked up ejecta material, and has elevated abundances of iron 
and intermediate mass elements (Si, S, Ar, Ca). ii) The donor 
star is still in the expanding phase: In this case the star should 
be underluminous and cooler, and may have the appearence of 
a later type star. Hi) The donor star has lost its entire envelope: 
In this case the donor's remnant should be a massive (~ 1 M ) , 
young CO white dwarf. 

5.3. Is the progenitor system of SN1604 typical for Type la 
supernovae? 

The model outlined here is specifically intended to explain many 
characteristics of Kepler's SNR. Nevertheless, given the many 
problems surrounding conventional models for Type la progen- 
itors (see the Introduction), it is worthwhile to discuss whether 
our model could be more widely applicable than just to this his- 
torical SN. 

The progenitor system in our model is a symbiotic binary, 
where the WD accretes mass from its AGB companion through 
non-conservative mass transfer. This progenitor system is not a 
conventional progenitor in the rich literature of Type la SNe. 
The reason is that RLOF of these binaries leads to unstable mass 
transfer, while wind accretion is characterized by low accretion 
efficiencies (< 10%), making it difficult to result in Type la 
events. However, these progenitor systems are currently still far 
from being understood. The Bondi-Hoyle model that is used to 
describe the wind ac cretion is ra ther simplified and not applica- 
ble to slow winds ( |Edgar|2004| ), while the much more efficient 
wind-RLOF mass transfer has not yet been studied in terms of 
population synthesis. Observations already seem to reveal sym- 
biotic systems with AGB donor stars, which, despite their large 
separation, show evidence of direct mass transfer and ^substan- 
tial hydrogen burning at the WD's surface (Mira AB : |Karovska 



etal.| ( |2004| ); V 407 Cyg: Mikolajewska 2010, private communi- 
cation). Furthermore, an non-negligible number of SNRs, iden- 
tified as Type la, display spectroscopical evidence of strong in- 
teraction betwe en the ejecta and dense CS M such as SN 2002ic 
(Ha muy et al.||2003]) and SN 20 05gj ( [Aldering et aT][2006l) and 
SN 20 06X flPatat et al.||2007) . In addition, |Borkowski et aL] 
( 2006 ) showed that the highly centered concentrated Fe emis- 
sion of the two Type la SNe in LMC DEM L238 and DEM L249 
can best be explained assuming that the ejecta interacts with a 
non-uniform CSM. All these observations point towards a SD 
scenario, where the circumstellar matter was formed by the slow 
wind of a giant, or an AGB donor star. Thus, non-conservative 
mass transfer of symbiotic binaries should be considered as a 
viable path towards Type la SNe. 

The generally favored SD model of a Roche-lobe overflow- 
ing donor star and an accreting, nuclear burning WD has two 
major problems with observations of Type la SNRs. First, the 
presence of neutral hydrogen around young Type la SNRs is at 
odds with the pre sence of large (~ 30 pc ) HII regions around 
supersoft sources ( [Ghavamian et al.|2003 ). Second, there seems 
to be a lack of t he observed X-ray emission in the elliptical and 
spiral galaxies ( [Gilfanov & Bogdan| ([20T0j>; |Di Stefano| ( [20T0| , 
respectively), which is expected from the steady hydrogen fusion 
on the WD surfaces (however see Hac hisu et al.|2010| ). The sym- 
biotic binary model for Kepler does not present either of these 
problems. 

The accreting WD in our model is enshrouded by the dense 
wind from the donor. The column density towards the WD in 
such a case is given by: 
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with R in corresponding to the separation of donor star and WD. 



Even for a mass-loss rate of M w = 3 x 10 6 M yr * we still 



have A^h = 10 cm -2 , sufficient to attenuate the U V/X-ray flux 
by a factor 10 4 (based on the absorption model of | Wilms et al.| 
2000 ). This does not only explain the lack of extended HII re- 
gions around young, nearby Type la SNRs, but it also suggests 
that the progenitors of SNe la may not necessarily be detectable 
in X-rays. 

Note that the mass-loss rates of the models in the range of 
10" 6 - 10" 5 M Q yr" 1 are lower than or near the upper limits in- 
ferred by the lack of narrow H and He emission lines and X- 
ray emission from SNe la ( [Mattila et al.||2005| |Hughes et al 



2007, respectively), but much higher than the strictest upper 
limit to mass loss rates set by radio o bservations, which g ives 
M ~ 3 x 10" 8 (u w /\0 km s" 1 ) M^r 1 ( [Panagia etal.|2006 ). 

The last limit poses a challenge to our model. However, even 
if the lack of radio flux in Type la SNRs puts clear constraints to 
the density of the ambient medium, the numerical lower limit es- 
timated by Panagi a et al. (2006 ) is debatable as it relies on semi- 
empirical models based on the radio emissivity of SNe Ib/c. In 
additio n, the solution suggested by |Wood-Vasey & Sokoloski| 
(2006), where successive nova explosions, just before the final 
Type la event, form an evacuated region around the explosion 
center, is also possibly applicable for the case of Kepler's SNR. 

An extr a challenge to our mo del is related t o the size of the 
donor star. Ha y den et al.| ( |2010| ) and |Tucker| ( |2011| ) analyzed 
the early-time lightcurves of several hundreds SNIa in the op- 
tical and UV band and they found no signatures from the blast 
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wave interaction with the donor star, predicted by Kasen ( [2010 ). 
This result constrains evolved stars with large radii as possi- 
ble candidates of WD companions in the SD Type la regime. 
Nevertheless, two reasons can possibly exclude our suggested 
model from this category. Firstly, even if the radius of a 4 - 5 
M© AGB star is large, a non-conservative mass transfer demands 
a wide binary system. Thus, the expected fraction of SNIa, which 
shows the signature of the collision at their lightcurve are much 
smaller (by a factor of ten) than the esti mated one for the Roche- 
lobe overflowing systems (see eq. 9 Kasen 2010). Secondly, 
most of the enhanced optical/UV emission caused by the in- 
teraction of the ejecta with the donor star will in our model be 
absorbed by the dense AGB wind. However, more quantitative 
conclusions regarding the light-curve for wind accreting Type 
la progenitors should be based on detailed hydrodynamical and 
radiative modeling of the immediate aft ermath of the Type la ex- 
plosion, similarly to what Kasen ( 2010 ) has done for Roche-lobe 
overflowing binary systems. 

A second implication of our proposed progenitor model is 
that the high column densities suggest significant absorption dur- 
ing the first few days following the explosion, when the super- 
nova blast is carving itself out of the denser parts of the wind. 
Assuming again M w = 3 x 10" 6 M Q yr" 1 , u w = 20 kms -1 , and a 
supernova shock velocity of 25,000 km s" 1 , then the optical light 
passes through a column of at least Nu ~ 10 22 cm -2 one day after 
the explosion, 4 x 10 21 cm -2 after four days, and 2 x 10 21 cm -2 
around 8 days after explosion, corresponding roughly with the 
time around maximum luminosity. Using the Galactic value for 
converting Nu into A v ( [Predehl & Sch mitt 1995 ), around peak 
luminosity the absorption is expected to be roughly Ay « 1. 
Therefore, the supernova should have had an appreciable intrin- 
sic reddening the first few days, if it originates from a symbiotic 



system. I nterestingly, early light curve data o f SN1994D (Patat 
|etal.|1996| ) and SN1996X ( [Salvo et al.|200T> , indeed, show that 
before maximum light the supernova initially becomes bluer and 
after maximum light becomes redder again. This might be ex- 
pected if the blue atmosphere is initially reddened by the wind 
bubble, and, once it has emerged out of the wind bubble, it be- 
comes redder again, but then caused by the adiabatic cooling of 
the ejecta. In this light it is also interesting that, when SN1604 
was first observed, its color was described to be "like Mars", 
and only on day eight it was described as "like Jupiter" (Baade 
|1943| ). Note that both Jupiter and Mars were in conjunction at 
the time, and SN 1604 appeared very nearby these planets, so at- 
mospheric extinction was very similar for all three objects. This 
conjunction was in fact the reason that SN 1604 was detected so 
early on ( [Green & Stephenson|2003| ). 

We like to end this section by pointing out that the one-sided 
shell, formed close to the explosion center of Type la SNe, could 
be rather common. The properties of the shell depend on the de- 
tails of the wind parameters, the systemic velocity of the progen- 
itor system and the local density of the ISM. Thus, even if the 
systemic velocity adopted here for the case of Kepler's SNR is 
quite extreme for low mass binary systems, similarly sized shells 
can be formed considering a progenitor system that has a much 
lower velocity but moves through a denser ISM. Interestingly, 
several Type la SNRs both show evidence for an axis of sym- 
metry, and seem to encounter a density enhancement in one di- 
rection. Two examples are Tycho's S NR (SN1572), which has a 
lower expansion velocity in the east ( [Katsuda et al.||2010 ), and 
SN 1006, which is very symmetric in the northeast/southwest di- 
rection, but has a higher density in the northwest. 

We offer these examples as tentative evidence that Kepler 
may not be as unique as we made it out to be in the Introduction. 



Rather, the one-sided shell in Kepler may just be more apparent 
than in other SNRs. 



6. Conclusions 

In this paper we have presented evidence that the observational 
characteristics of the Type la SNR Kepler/SN1604 are best ex- 
plained by a moving progenitor system, in which the donor 
star forms a wind blown bubble, while at the same time the 
WD accretes part of the wind, thereby reaching 1.38 M within 
~ 0.4 - 0.7 Myr. Due to the systemic velocity the wind bubble 
is surrounded by an asymmetric shell, reaching its highest den- 
sities in the direction of the motion of the progenitor. In order to 
explain the nitrogen abundance in the north of Kepler, the donor 
star probably had a main sequence mass of 4-5M©, which are 
known to have nitrogen rich envelopes in the AGB phase. 

We used hydrodynamic simulations to show that the wind 
properties of the system can explain the observational charac- 
teristics of Kepler's SNR, namely a one-sided shell with which 
the SNR is interacting, and a slower expansion velocity in that 
region. 

Our simulations show that in the direction of the progeni- 
tor's velocity the SNR blast wave may just have completely pen- 
etrated the wind blown shell, or still be inside it. Based on the 
possible presence of X-ray synchrotron radiation in the north of 
Kepler, and the lack of it in adjacent regions, and the absence of 
nitrogen rich knots in that direction, we argue that in Kepler part 
of the blast wave may indeed have penetrated all the way through 
the shell. On the other hand, the presence of B aimer dominated 
shocks indicate that at least part of the blast wave must still be 
within the cool, neutral, wind shell. 

Our results also show that the distance to Kepler's SNR is 
likely to be around 6 kpc, and not 4 kpc ( Sankrit et al. 2005 ). This 
larger distance estimate is in agreement w ith recent values based 
on the absence of gamma-ray emission ( [Berezhko et al.|[2006} 
Aharoni an et al. 2008| ), and the kinematics of Kepler's SNR 
( Vink 2008 ). However, a distance of 4 kpc can be reconciled with 



the simulations if the explosion energy was 2 x 10 5U erg. This 
would make SN 1604 a subenergetic explosion, which seems 
unlikely given its historical light curve and the copious amount 
of Fe present in the SNR. According to our scenario, the donor 
star should still be present within the SNR. It is likely an evolved 
star with most of its envelope mass stripped. 

We have argued that our model may be more broadly appli- 
cable to Type la supernovae and their progenitors. The presence 
of a dense wind close to the supernova should result in consid- 
erable intrinsic reddening of the supernova during the first few 
days after the explosion. For the pre-supernova appearance of 
the system, the large absorption column blocks out most of the 
X-rays from the accreting WD. These systems would, therefore, 
be classified as symbiotic binaries, and not as supersoft sources. 
This also explains the presence of neutral hydrogen near young 
Type la SNRs, which cannot be present around supersoft sources 
due to their high X-ray/UV flux ( Ghavami an et al.|2003| ). 

Another consequence of the model, if it is more widely ap- 
plicable to Type la progenitors, is that many SNRs should be 
interacting with one sided shells. For Tycho's SNR this could 
explain the low expansion rate in the east, and even, at least par- 
tially, the proximity of the ejecta close to the forward shock, an 
observational characteristic usua lly attributed to the presence of 
cosmic rays ( [Warren et al.|2005 ). 

Kepler's SNR requires a rather narrowly defined model. So 
major question remains whether SN 1604 was a very special 
event with an unusuale progenitor history, or on the other should 
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be considered to be a key object for a better understanding of 
Type la supernovae in general. 
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